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The material presented here is supplementary to the article entitled: "Density-functional investi- 
gation of the rhombohedral to simple cubic phase transition of arsenic" ; it deals with the convergence 
issues involved in studying a semi-metal to metal structural phase transition such as the A7 sc 
transition of arsenic. 



The occurrence of the A7 sc phase transition of ar- 
senic is identified most clearly by the behavior of the 
nearest and next-nearest neighbor distances. In this 
study we investigate closely the convergence of these dis- 
tances with respect to fc-point grid size and smearing over 
the pressure ranges of 12-24 GPa and 17-30 GPa for the 
LDA and GGA-PBE cases, respectively. In the region 
of the transition, the energy differences between the two 
structures are extremely small. If for a given pressure 
the potential energy surface is noisy and rather flat, then 
it is possible that our system relaxes to a structure that 
corresponds to a local minimum, rather than the global 
minimum of our energy surface. We have found that 
to be able to achieve good results we must resolve our 
potential energy surface as much as possible. We have 
investigated Gaussian smearing^ and Mcthfessel-Paxton 
smearing)^ but we find that we obtain much cleaner and 
clearer results by using cold-smearing^ and therefore have 
done so throughout this study (a review of these different 
smearing techniques can be found in Rcf. 0). 

In these numerical calculations, the integrations that 
are to be performed over the Brillouin zone are dis- 
cretized by way of the fc-point grid. In the case of an insu- 
lator, no smearing is required. Smearing is only required 
in the event that one or more bands cross the Fermi level, 
in other words if the material possesses a Fermi surface 
(if it is a metal or a semi- metal). The effect of smearing 
is to blur the details of the Fermi surface (by imposing an 
artificial temperature on the electronic system); smear- 
ing must be used to ensure that the calculation is stable, 
but can also be used to hasten the convergence of the 
calculation. If our material does possess a Fermi surface, 
but we use an infinitely dense fc-point grid, no smearing 
is needed as no details of the Fermi surface will be lost; 
the choice of fc-point grid and smearing are coupled to 
each other. We would like to minimize the number of 
fc-points used in the calculation, but must choose a grid 
dense enough to properly sample the Brillouin zone as 
well as to pick out the features of the Fermi surface. If 
we use a smaller (less dense) fc-point grid, we must use 
more smearing, but the exact amount of smearing to use 
is not obvious. If we use a smearing that is too high, 
details of the Fermi surface may be washed out possibly 
to the extent of affecting the quantities that interest us. 



It is preferable to keep the smearing low to obtain re- 
sults of the highest accuracy possible, but there is a limit 
to how low the smearing can be set for a given fc-point 
grid before there is a risk that the calculation becomes 
unstable. 

In fact, to choose the smearing properly for a partic- 
ular fc-point grid would require knowledge of the three- 
dimensional band-structure of the material throughout 
the entire Brillouin zone, since the smearing should be de- 
termined by the steepest gradients of the bands crossing 
the Fermi level. (Although we cannot avoid it, it is not 
actually appropriate to apply a unique amount of smear- 
ing to the entire band-structure of a material undergoing 
a semi-metal to metal phase transition; ideally we would 
use some sort of adaptive smearing technique to minimize 
the washing-out of the features of the Fermi surface.) 
It is not feasible to investigate three-dimensional band- 
structures, so we must carefully converge our calculations 
with respect to both fc-point grid size and smearing. Ini- 
tially, we did some cursory convergence tests to determine 
suitable choices for the cut-off energy as well as for the 
density of the fine-grid, respectively to ensure an appro- 
priate basis set size and to recover sufficient detail of the 
charge density within the atomic cores; these tests were 
performed on the uncompressed system. The results of 
our geometry optimizations, however, must be subjected 
to more rigorous convergence testing; each calculation 
must be repeated in its entirety for each combination of 
fc-point grid and smearing. The properties that interest 
us must be converged in this way before anything can 
reliably be said about them. Although convergence tests 
for fc-point grid size are standard, in electronic structure 
calculations that involve smearing it is very rare to find 
that testing has been undertaken to determine how the 
smearing has affected the results; yet applying a casual 
approach to the choice of the amount of smearing used 
can lead to gross inaccuracies (see Ref. i). Thus, it is 
essential that tests for convergence with respect to both 
fc-point grid and smearing be performed. Furthermore, 
it must be stressed that such tests should not merely be 
carried out by looking at what happens to the total en- 
ergy of the system; rather, these tests must always be 
carried out on the quantities specific to the study itself. 

As we mentioned above, we would like to perform an 
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in-depth investigation of the A7 — ^ sc phase transition 
and so we must test the convergence of our quantities of 
interest, the nearest and next-nearest neighbor distances, 
with respect to fc-point grid and smearing in the region 
of the transition. Fig. [T] displays the results of these con- 
vergence tests for the case of the LDA. Consider first the 
top panel of this figure, which displays the behavior of 
the nearest and next-nearest neighbor distances for dif- 
ferent grid sizes using a cold-smearing of 0.1 cV. We see 
that for this value of the smearing, none of the results 
for the different grid sizes have converged in the imme- 
diate vicinity of the transition, which appears to occur 
somewhere in the range of 19-21 GPa. Even our results 
for a 50x50x50 grid have not quite converged to those 
for a 66x66x66 grid. Note also how different our results 
at the transition arc for the 24x24x24, 25x25x25 and 
26x26x26 grids at this smearing. This is down to the 
particular sampling of the Brillouin zone for each specific 
grid and tells us that we arc not using enough fc-points 
to properly sample the Brillouin zone for this smearing 
and at these pressures. Outside the region of the phase 
transition, the calculations are less sensitive to the den- 
sity of the fc-point grid; below about 14 GPa and above 
about 30 GPa we see that our results are well converged 
at this smearing for all grid sizes. Hence we can be confi- 
dent about our choice of grid and smearing (33x33x33, 
0.1 eV) at pressures away from this particular transition. 

If we decide that we do not require the finest features of 
the Fermi surface to obtain a meaningful result, then we 
can increase the smearing, as we have done for the middle 
panel of Fig. [I] for which a cold-smearing of 0.2 eV was 
used. We see that at this smearing, the results have con- 
verged for the three most dense grids used, 33x33x33, 
50x50x50 and 66x66x66, and the transition pressure 
seems to have shifted slightly higher appearing now to 
occur at 20-21 GPa. 

In the lowest panel of Fig. [Tl we have increased the 
smearing high enough to ensure that the results of all 
of our grids converge. Having shown, for a smearing of 
0.5 eV, the convergence of the calculations for the least 
dense grids used, we do not repeat the calculations for 
the most dense grids of 50x50x50 and 66x66x66. If we 
compare the top and bottom panels, we can see quite 
clearly that as the smearing is increased, the transition 
shifts to slightly higher pressures. 

The components that make up the nearest and next- 
nearest neighbor distances are the cell angle a, the atomic 
positional parameter z, and the lattice parameter a. We 



can therefore examine the convergence properties across 
the transition of these constituent quantities as well. In 
the case of the LDA, the results of our convergence test- 
ing on a, z, and a are presented in Fig. [21 Fig. [3] and 
Fig. m respectively. 

We observe similar difficulty in the convergence of the 
A7 sc phase transition in the case of the GGA-PBE. 
Fig. O shows the results in terms of the nearest and next- 
nearest neighbor distances of our GGA-PBE calculations 
for each of the same grids employed for the LDA, but for 
a cold-smearing of 0.1 eV only. Considering the data 
resulting from the use of the most dense fc-point grids, 
we can conclude from this figure that in this case the 
phase transition must happen over the pressure interval 
of 27-29 GPa. The corresponding results of the conver- 
gence testing on a, z, and a in this case are revealed in 
Fig. [6l Fig. [7] and Fig. [U respectively. 

As we said above, for a cold-smcaring of 0.1 eV even 
calculations using a 50 x 50 x 50 fc-point grid have not con- 
verged to those using a 66 x 66 x 66 grid in the vicinity of 
the transition. We can conclude just by observing the 
top panel of Fig. [T] that in the vicinity of the phase tran- 
sition, these calculations are extremely sensitive to the 
details of the Fermi surface. The Fermi surface of arsenic 
is extremely complex. From work carried out by Lin and 
Falicov^ we estimate that a fc-point grid at least twice as 
dense as our most dense grid used (66x66x66) would be 
required to resolve all of the features of the Fermi surface 
of arsenic at ambient pressures. 

The most dense grid employed by Da Silva, et. al^ in 
their study of arsenic was 13x13x13, and in the case of 
Haussermann, et. al^ it was 17x17x17. Durandurdu^ 
uses only the gamma point for a unit cell containing 250 
atoms, roughly corresponding to using a 5x5x5 fc-point 
grid for for a two-atom unit cell. Our investigations re- 
veal that these calculations could not have been con- 
verged, and that any agreement with experiment that 
may have been observed in these cases was merely fortu- 
itous. 

We conclude that it is surprisingly difficult to converge 
the calculations for this semi-metal to metal phase transi- 
tion using reasonably sized fc-point grids and smearings. 
It is of course possible to converge the calculations us- 
ing less dense grids if very large smearings are used but 
it would be at the expense of accuracy in the resulting 
transition pressures. To ensure accuracy when studying 
a pressure-induced semi-metal to metal phase transition, 
dense fc-point grids are essential. 
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FIG. 1: Nearest neighbor (solid lines) and next-nearest neighbor (dotted lines) distances as a function of pressure for the A7 — > sc 
transition of arsenic using the LDA. Grids used are: 24x24x24 (black), 25x25x25 (red), 26x26x26 (orange), 33x33x33 (cyan), 
50x50x50 (blue) and 66x66x66 (green). The values of cold-smearing used are: 0.1 eV (top), 0.2 eV (middle) and 0.5 eV 
(bottom). 
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FIG. 2: Cell angle a as a function of pressure for the A7 — > sc transition of arsenic using the LDA. Grids used are: 24x24x24 
(black), 25x25x25 (red), 26x26x26 (orange), 33x33x33 (cyan), 50x50x50 (blue) and 66x66x66 (green). The values of 
cold-smearing used are: 0.1 eV (top), 0.2 eV (middle) and 0.5 eV (bottom). 
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FIG. 3: Atomic positional parameter 2 as a function of pressure for the A7 — > sc transition of arsenic using the LDA. Grids 
used are: 24x24x24 (black), 25x25x25 (red), 26x26x26 (orange), 33x33x33 (cyan), 50x50x50 (blue) and 66x66x66 (green). 
The values of cold-smearing used are: 0.1 eV (top), 0.2 eV (middle) and 0.5 eV (bottom). 
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FIG. 4: Lattice parameter a as a function of pressure for the A7 — > sc transition of arsenic using the LDA. Grids used are: 
24x24x24 (black), 25x25x25 (red), 26x26x26 (orange), 33x33x33 (cyan), 50x50x50 (blue) and 66x66x66 (green). The 
values of cold-smearing used are: 0.1 eV (top), 0.2 eV (middle) and 0.5 eV (bottom). 
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FIG. 5: Nearest neighbor (solid lines) and next-nearest neighbor (dotted lines) distances as a function of pressure for the 
A7 sc transition of arsenic using the GGA-PBE. Grids used are: 24x24x24 (black), 25x25x25 (red), 26x26x26 (orange), 
33x33x33 (cyan), 50x50x50 (blue) and 66x66x66 (green), each with a cold-smearing of 0.1 eV. 
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FIG. 6: Cell angle o? as a function of pressure for the A7 — > sc transition of arsenic using the GGA-PBE. Grids used are: 
24x24x24 (black), 25x25x25 (red), 26x26x26 (orange), 33x33x33 (cyan), 50x50x50 (blue) and 66x66x66 (green), each 
with a cold-smearing of 0.1 eV. 
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FIG. 7: Atomic positional parameter z as a function of pressure for the A7 sc transition of arsenic using the GGA-PBE. 
Grids used are: 24x24x24 (black), 25x25x25 (red), 26x26x26 (orange), 33x33x33 (cyan), 50x50x50 (blue) and 66x66x66 
(green), each with a cold-smearing of 0.1 eV. 
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FIG. 8: Lattice parameter a as a function of pressure for the A7 — > sc transition of arsenic using the GGA-PBE. Grids used 
are: 24x24x24 (black), 25x25x25 (red), 26x26x26 (orange), 33x33x33 (cyan), 50x50x50 (blue) and 66x66x66 (green), 
each with a cold-smearing of 0.1 eV. 



